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摘 X 本 文 基于 变 分 多 尺度 方法 提出 了 非 定常 Navier-Stokes 方 程 的 一 种 稳定 化 方法 ， 并 利用 一 个 与 实 
际 误差 等 价 的 后 验 误 差 估计 子 ， 结 合 自 适应 算法 ， 得 到 了 非 定 常 Navier-Stokes 方程 的 自 适应 变 
分 多 尺度 稳定 化 方法 。 该 方法 简单 、 直 观 ， 且 易于 程序 实现 。 数 值 算 例 验 让 了 该 方法 的 有 效 性 。 
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有 限 元 方法 是 Courant 于 1943 年 首先 提出 。 二 十 世纪 60 年 代 ， 有 限 元 方法 的 应 用 范围 扩展 
到 了 结构 力学 以 外 的 领域 ， 并 建立 了 严格 的 数学 理论 基础 ， 创 立 了 相应 的 商业 软件 包 。 直 到 现 
在 ， 有 限 元 方法 已 经 形成 了 比较 完备 的 理论 体系 ， 成 为 工程 中 应 用 最 为 广泛 的 数值 方法 。 在 采 
用 有 限 元 方法 求解 具体 问题 时 ， 其 网 格 剖 分 的 如 何 对 计算 结果 的 精度 起 着 决定 性 的 作用 。 正 因 
为 如 此 ， 近 年 来 出 现 了 基于 有 效 误差 估计 的 自 适应 有 限 元 方法 。 

自 适应 有 限 元 方法 的 思想 最 早出 现在 1978 年 ， 美 国 数学 家 Babuska 等 针对 线性 椭圆 问题 在 
文献 [1] 中 提出 了 一 种 数学 上 严格 的 基于 残 值 的 误差 估计 技术 ， 给 出 了 误差 的 界限 ， 并 用 局 部 
误差 指标 定义 了 最 优 网 格 。80 年 代 以 后 ， 自 适应 有 限 元 方法 方面 的 研究 成 果 很 多 2a ， 发 展 极 
为 迅速 。 目 前 ， 自 适应 有 限 元 方法 有 h- 加 密 方法 、p- 改 进 方法 、 移 动 结 点 法 、h-p 法 、r-p 法 
等 。 

本 文 基 于 文献 台中 的 变 分 多 尺度 方法 的 基础 上 ， 考 虑 由 非 泡 函数 扩充 的 分 片 线性 有 限 元 空 
间 ， 取 线性 有 限 元 为 粗 尺度 ， 非 泡 函 数 为 细 尺 度 ， 从 而 提出 了 Navier-Stoke 方 程 的 一 种 稳定 化 
方法 。 此 外 ， 结 合 文献 [5] 中 Berrone 提 出 的 白 适应 算法 ， 我 们 给 出 了 两 个 后 验 误差 估计 子 ， 并 
理论 证 明了 其 与 实际 误差 的 等 价 性 。 

本 文 内 容 安排 如 下 : 第 二 部 分 概括 介绍 了 非 定 常 Navier-Stokes 方 程 的 标准 有 限 元 解法 ; 第 
三 部 分 介绍 了 非 定常 Navier-Stokes 方 程 的 变 分 多 尺度 方法 及 其 稳定 性 ; 第 四 部 分 与 第 五 部 分 
介绍 了 自 适应 算法 以 及 后 验 误差 估计 子 与 实际 误差 的 等 价 关 系 ; 第 六 部 分 给 出 了 经 典 的 方 腔 流 
实验 算 例 ， 从 而 说 明了 本 文 方法 的 有 效 性 。 
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* 基 金 项 目 : 国家 自然 科学 基金 (10871156); 西安 交通 大 学 交叉 学 科 项 目 基金 (2009xjtujc30). 
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2 AEH Navier-Stokes 方程 的 标准 有 限 元 解法 


考虑 非 定常 Navier-Stokes 方 程 
ди — „Ди + (u: V)u + Vp = f, (x,t) € Q x (0, T), 
У.и=0, (x,t) € Q x (0, T), 
и = 0, (2,1) є AQ x (0, T'), 
u(-,0) = 0, z€ Q, 
这 里 Btu = 95 QER rh ЙН s SWB Rh, Hi OQ E Lipschtiz E & f , fe £ 
Av > 0, WHT > 0, SERR u = Ga(z,t),u2(z,t)), RAp = p(z,t), HF = f(z,t) € 
C9(0,T; [L7(Q)]?). 
首先 ， 我 们 定义 


Hi(Q) = {v € H1(Q) : vlan = 0}, (9) = (ae £49): f 4-0) 


(1) 


以 及 V = HiQ), Q = LO. WC), | - PRA Si [LO], d = 1,2 Ef P2 Bz 
MADR, |: h 与 | :|z 分 别 为 空间 V 上 的 HH! 半 范 与 HH?2 半 范 。 为 了 得 到 原 问 题 (1) 的 
标准 有 限 元 解法 ， 我 们 把 时 间 (0,T) 分 成 N 个 时 间 段 ， 每 个 时 间 段 记 为 (t"-1,t")， 其 步 长 
Эт" =e, sn<N, Rih0= 0 < <... ct" T. ABR D, Wu = 
u(t) р" = p(. tn), fm = f(. tn). RIEV x V xV 上 定义 三 线性 项 


шош) = 2 ((u- Vo, ш) — 2 ((u* Vows»), 
AA OA F| EE {ЕД 
b(u,v,w) = —b(u,w,v), Vu,v,w € V, 
би, v, w)] + |b(o, и, w)| + lbw, и, v)] < сш. о, 
VueV, veD(A) = [H2(Q)2 n V, we [22(9))2, 
那么 原 问题 (1) 的 变 分 形式 为 ， 求 (wu",p") eV x Q, ЕЖ (u, q) € V x Q # 
Во([и", р"), [v,g]) = (vVu^, Vv) — (p^, V v) + (q, V - u”) 
= (f",v) — (@iu|i= v) — Би", u^, v), 
В Во([и", p^], [v, q]) WE: 存在 一 常数 8B， 对 任意 的 (wu",p")eEVxQ@ 有 


O< (uh lp ll) < sp Brod) 


(2) 
(v,q)eV xQ lv]: + llall 


FANT BC (TP), 5 n 4 BJ [н] Ez P< 2 Q BJ ba EAR ЖИЛ, SEA A RA ЛЕМ $É £ 1 
Aeg, ЖШ F BJ 48 BR лол [a] 


Ve = {op e V NOMP :vlr e [PK], v K e Tr}, 


Qh= (qh ERNA) : ак e PK), VK e TR}, 
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其 中 POCO) 是 次 数 i > 1 的 多 项 式 构成 的 空间 ， 记 
hy = diam{K : K € Ty, h=max(hx), he ={|E|: E € ef}. 
于 是 ， 我 们 得 到 了 原 问 题 (1) 的 全 离散 变 分 形式 ， 对 任意 的 n= 1,…,N, (u, qh.) € Vp x Q, 
Ж (ug ph) € V x Q 使 得 
Bi( [wR, Ph]; [vn, qh] ) + b(uy, upt, Un) + b(un 1, ti Un) 
= (f",Un) + b(ug up’, vn), (3) 
这 里 
u? — utl 
Bi (luk, vn), [of an] ) = (At, vn ) + (Vuk, Von) — (PR, V n) + (ans V - ub). 
众所周知 ， 当 我 们 采用 和 低 阶 等 次 有 限 元 对 Р1Р1 (或 Q1Q1) 求 解 问题 (3) 时 ， 得 到 的 解 是 不 


稳定 的 。 许 多 学 者 相继 提出 了 多 种 稳定 化 方法 813。 下 面 我 们 就 考虑 用 变 分 多 尺度 方法 来 求解 
问题 (1)， 以 达到 稳定 化 的 目的 ， 使 其 得 到 稳定 解 。 


3 ”求解 非 定常 Navier-Stokes 方程 的 变 分 多 尺度 方法 


利用 原 问 题 的 离散 变 分 形式 求解 稳定 的 数值 解 时 ， 变 分 多 尺度 方法 是 一 种 常见 的 稳定 化 方 
法 。 变 分 多 尺度 方法 一 般 可 以 通过 两 个 步骤 来 完成 。 

步骤 1 ”把 原 问题 的 数值 解 分 成 两 个 部 分 ， 粗 尺度 部 分 与 细 尺 度 部 分 ， 相 应 的 检验 函数 也 分 
成 了 两 个 部 分 ， 然 后 把 检验 函数 的 两 个 尺度 部 分 分 别 代入 到 原 问题 的 离散 变 分 形式 中 ， 这 样 原 
问题 的 离散 变 分 问题 化 成 了 两 个 问题 ， 粗 尺度 问题 与 细 尺 度 问题 ; 

步骤 2 从 步骤 1 中 的 细 尺 度 问 题 中 求 得 数值 解 的 细 尺 度 部 分 与 粗 尺度 部 分 的 关系 ， 然 后 把 
这 一 关系 代入 到 粗 尺 度 问题 中 ， 从 而 使 问题 得 到 了 解决 。 

类 似 于 文献 加 中 的 方法 ， 假 设 细 尺 度 部 分 关于 时 间 t 是 分 片 常数 ， 我 们 可 以 得 到 原 问 
题 (1) 的 变 分 多 尺度 方法 ， 对 任意 的 (ww gn) € V? x Ой, R (un, pp) € V? x Q? 使 得 


ип — ug! - - 
( 8—28— on ) + Bo([ug, pR], [vns an] ) + (up ug, vn) + bluR t, uR, vn) 


hof n up! а n-1 ,n-1 
= У (7 + zn Va) + (f ‚®һ) + b(up UR o Un); (4) 


Jtrh 


h% pu 2 7 
+ У c + Уру, Van) + (ans V - up) 


h 
+ > om [одлар в, [vO, un] E )E 
EceR 
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3.1 稳定 性 分 析 


定义 如 下 范 数 
lluk П vius + Y 2 ell! (рда) |2, 12 = Y. КЕР x, 
ЕС! KET 
我 们 得 到 下 面 的 结论 


定理 1 BV, Qh) eMail, MER (uh ph) € Vp x Оһ, (u",p"), („4 eV x Q, fé 
在 与 无 关 的 常数 c 和 6， 使 得 当 h < cr* 时 ， 有 下 列 式 子 成 立 


Bo( [uf pi], [uf pp] = wh IÈ ПРЕПА + 2. ME (SE ype) ， (5) 
| Ba (u^, р"), lo, oD)| < e(lu”h + р") (els + lal); (6) 
о< (иһ +1 1)< sp Zei peh aD (7 


(wonanevexan — lvnh + llaa || 


证 明 由 Bo se X RT A (5) 式 成 立 。 另 外 ， 我 们 有 局 部 迹 定 理 ， 存 在 与 疡 无 关 的 常数 c > 
0， 使 得 


1 
wl < e ( gz lol + halle). voe [I G0)", (8) 
由 Schwarz 不 等 式 知 
hg 
TE (panule, oro) )s < FE | owls |1 деје |, 
而 由 逆 不 等 式 与 (8) 式 知 


hg ni 112 oe 
951194") Ile 


I^ 


Top ^K Hip vu? 112. + hz |u Vu? п к) 


cvhg, n cvhg, n 
he = lu Le + тре ju фк < Sle k 


I^ 


所 以 


hp " " 
xm (Inu le, [идлаје) p < clu" hxloh xc 


ШТ h < ст", t Schwarz К |а| € clues 知 


hi. pu” hk уол 
— < n 
Mic (gg) < E E Lane рачка сит lalz, 
所 以 

"n n TL ul n hig u” n 
| Bo([u p ], [v, a) | x v(Vu ,Vv) — (p ‚„М.0) + (9, V- ug) 十 x = (+ Va), 

KET; 
< vju^llel + "Ы + lalek + > еа, + lip" lace) 


KETE 


I^ 


с(1ш" + [lp ll) (lela + llall), 
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ШП (6) 式 得 证 。 
最 后 ， 我 们 证 明 (7) 式 成 立 。 我 们 的 目的 是 寻找 一 对 (vy gs)， 使 得 


Вә ( [ив, рь], vn, qa] ) > є(|ир | + П) (Juni + Ilall ). 


首先 ， 我们 令 (vn, аһ) = (—w, 0) (и 待定)， 利 用 Cauchy-Schwarz 不 等 式 ， 有 


B (lun, ph], [-w,0]) = ~v(Vul, Vw) + (nf, V -w) — F e оди], [vs] j 
Eceg 
> -vho - У 28| д,ар) ell baw] [lg + GS V) 
ECer 
> (wii E 0) (эш + У r leali)’ 
ЕСєр ECe? 
+(р?, V w) 
> (PR, Vw) ll ug Ila (ш + V7 ae | Ponal lla) (9) 
ЕЄє 


由 (2) 式 知 ， 对 任意 的 成 e QR, FFE w e V 8 wl = |р, Н(У-ш,рр) > сіш р. F 
是 ， 由 上 式 及 (8) 式 与 (9) 式 ， 我 们 有 


Вә([иҗ,›рь],[—ш,0]) 2 —Vev || uË llla Jwh + сі): lip? 

—evy,! il vg (12 +e- ә), (10) 

Ж y 充分 小 (待定 )。 我 们 取 (vn, an) = (un 一 бш,рь), ó € R2 (待定 )， 由 (5) 式 与 (10) 式 知 
Ba ( [ug ph], [on, as] ) 


= Ba( [ Р, ph]; luk pp ) + óB ( [luk ph], [~w, 0]) 


IV 


Е hi, / un 
> (1 волу) It ME +e- vlone tote + Y^ PEE vy) 
KeTp 
_ n n h? hi ир п 
zov ~ bevy, *)jub|? + (e — т) + > рв к + у? З (Sven) an) 
KET? KET? 


另 一 方面 ， 由 Cauchy 不 等 式 知 


hk (UE ga hk hi 
ee ny2 _ “K n2 . 
8v (2 ; VPh jm = 16k) lus lc l6v Pili, к 


由 于 [ug]  cjug|,, ЖШ 
hi UR n C 2 hi 
UK 人 “h > oca — lpn? ,.. 
Sv (35, ver) 2 ie li 16r РЕК 


H (11) Rat, 只 须 0 < < c, 0<ó< giller Säl: 有 


Ba( [uj ppl, оа) > v(1— 50 — 


c n n 
5 luni? + б(с — у) 150, 
16v 
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我 们 取 z i 
e(v) = min { E (1 бои? — 55), 396-9) }, c(ó) = 1-4 ó, 


则 
B>( [un ph), vn. an] ) > с(и) lug] + lips? ). 
[бл] + lanl] < lu] + [|р + бш < jul Q + ó)||p# < с(8) jul + 1211), 
WG = Z, M (7) RME, 证 毕 


а ” 自 适 应 算法 


在 文献 [4-6] 中 ， 作 者 已 经 指出 数值 解 的 离散 误 其 主要 有 两 个 来 源 ， 一 个 是 空间 离散 产生 的 
误差， 一 个 足 时 间 离 散 产 生 的 误差 。 为 了 较 好 的 控制 这 两 个 误差 米 源 ， 我 们 采用 下 向 两 种 误差 
估计 子 ， 并 结合 白 适 应 算法 ， 这 样 使 我 们 的 问题 得 到 了 较 好 的 解决 。 

4.1 空间 自 适应 

从 前 面 的 分 析 过 程 来 看 ， 我 们 发 现 空间 离散 误 鞠 主要 米 源 于 细 人 尺度 部 分 ， 也 就 是 说 wj 是 主 
要 的 误 莽 米 源 。 芬 外， 考虑 到 不 可 压缩 性 ， 我 们 引进 下 面 空间 局 部 的 误 革 估计 子 
ий — up! j 

r^ K 


n 


hl us 
nk = | || fe у - 


1 
: ; 
+ 2:3 Sell би |5. + IE V uz lla | 


其 中 fp fn BA TTE BOB R, RIENE = vius? + [ру [2 29 T Е® ТЫН Et L48312 
适 的 网 格 剖 分 ， 我 们 采用 均衡 误 鞠 估计 子 的 方法 ， 即 要 求 在 每 个 单元 上 ， 都 有 下 奋 的 不 等 式 成 
КУЛ 


(1 = o)? 8t (E)? 


Nx < (nk)’, (12) 
a 2 ^2 
(m )2 < GERE. (13) 


ДП oi, 81 是 给 定 的 常数 ，NEk ERE УНУ А, Ha € (0,1]， 称 [Bi 为 空间 上 的 
容忍 误 产 。 当 菜 个 时 间 段 闫 定时 ， 在 每 个 单元 上 ， 我 们 要 求 (12) 式 与 (13) 式 同时 成 立 。 
E (7&)? 比 (12) 式 左 端 的 下 确 界 小 时 ， 我 们 需要 放松 网 格 ， 相 反 ， 若 (02)? ЊЕ (13) 式 右 端的 上 
确 界 大 时 ， 我 们 加 密 网 格 。 经 过 加 和 密 或 放松 网 格 后 ， 我 们 在 新 的 网 格 上 求解 问题 (4)， 如 此 反 
复 下 去 ， 直 到 在 每 个 单元 上 (12) 式 与 (13) 式 都 成 并 。 

Ж 当 记 减 小 时 ， 网 格 剖 分 的 单元 总 数 增 加 ， 同 时 空间 上 上 解 得 精度 也 会 提高 。al 是 用 
KSB (nk) ME. FAR, ay = 工时 ，( 吸 ) 无 下 界 ， 网 格 剖 分 不 会 出 现 放 松 的 情 
№. oi, Bi 取 值 不 同 会 导致 计算 结果 的 人 不同 ， 在 这 里 我 们 个 考 虐 这 种 影响 情况 。 

4.2 时 间 自 适应 

根据 文献 [5] 中 的 分 析 结 果 ， 我 们 采用 以 下 时 间 局 部 的 误 美 估 计 子 


ту = yvr” [ug — REP 
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同样 ， 引 进 变量 on” = ТТЕ, ESB RER FIRER 

(1 - оз) < TÈ < (1+ оз), (14) 
X rh оз, ВЛЕН, Паз є (0,1, 8 Wy 22 Tal LAA LR OE. AE TT КОЙ 
足 (14) 式 ， 我 们 将 进行 下 一 个 时 间 步 的 计算 ， 理 则 ， 根 据 (14) 式 ， 我 们 需要 增加 或 缩小 时 间 步 
长 ， 直 到 在 该 时 间 段 上 (14) 式 成 立 。 我 们 同样 也 不 讨论 ao, В, 取 不 同 的 值 对 计算 结果 的 影响 情 
Who 


5 ”整体 误差 估计 子 与 实际 误差 的 等 价 性 证 明 


当 1 <n < N FF, 我们 定义 下 列 变量 


e= Y -ulk e = > 10" -rlz 
кєт KeTr 


m= | GR. e3 = \т"(ет+еў), 
KET? 


并 定义 整体 误差 估计 子 Н" BERRA e" 
н" = |Y (т"(ү")? 9»), e= |Y em. 
m=1 m=1 


我 们 有 下 面 的 结论 成 立 。 
定理 2 (ur, pe) є VP x Q 是 问题 (4) 的 解 ， 则 


n 
H^ «el Ут" у h(I = fe + un lun к 
m=i 


кєт 


二 п 
up т. up. 2 n m 一 ue 
+ | ú Әг") 十 (en)2 + ти" — y" 8} : (15) 


Tn 


n m m-1 2 
= up —u 
e он") + уст” у [|P - beulen |. + ARIE” et] 
m=1 кєт” 
z m m2,,,m|2 -2 m2,,,m—1,2 —1,4 1/2 
+ т" [fumes + a? (lulu! + fume } . 


m=1 


(16) 
证 明 令 


Ei; = | Ак, Ак,, ак = У up’, 


hg 
шк = У sr Фе[идһит]е 
ЕСӘК 24v 


m-—1 


h? m. 
- a(a){ SE (f -vop h) + b»: 25 (và, u]eer }, 
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ІХ, Ак,, Ак,, Ак, 是 单元 五 的 重心 坐标 ， 函 数 a(z) 定义 如 下 


1，2ZE 五 ， 
a(x) = 
| 0, z€ ƏK 
从 而 我 们 得 到 
m-1 


Е [vnut | =—-К{(ү{т_ хут иһ up 
K = si 0nun]g, wklk = gp V^ — Vph = K 
以 及 
h? um — ym! 2 hg 2 
(ng)? = Ee- veg | + $^ sZlló,uple|| + IS - url 
ECOK 
ир = uy m m m 
= (A= wx) | — (Ver, Vuk)k + (ez, V-wk)k — (UV -єү,чк)к 
+(f™ — fi’, wK)K — (дш ҥ=ҥт,шк)к — (u^ - V)u™ ,wk)k. 
rh ЬЕ! 
—((u™-V)u™, wx) < clu", x|u"l, x wx lie; 
所 以 
(ng < e| (16 + 21) (ек + ЇчкЇк) 
m m-—1 
m m u, 一 也 
+ (urn = flc | A - aule]. Im кїн” c) ll]. 
Y R sel (ep) er" Y^ AR (IS - Fk 
кєт," KeTp 
-1 
m2 m2 Up i Up. =. | 
+ u^ как + [z дри X 
另外 
(ng = rm 加 一 如- 性 <c[(eg)2 + lu" — u], 
所 以 


n 


(H" < ef Ут" Y aI- A + Le lB к 


т=1 кєт" 


um — ym-i 2 Ш 
JER аар) er merce] 
m=i 


Tm 


BU (15) лу. 
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为 了 证 明 (16) 式 成 立 ， 我 们 有 
Bo([u^, p^], [v,g]) = (vVu^, Vv) — (p^, V v) + (q, V ш") 
= (f^,v) — (Oli v) — ((u" - V)u^, v), 
Arh (4) XA 


B. (lun, ph], [оһ, qa ]) = (Уик, Von) — (ph. V - vn) + (qn; У ип) 


hie ү Uh h n 
E “к лш 
h 


i 
< 
a 
e 
= 
+ 
МА 
— 


—b(uz-!, uf, vx) 一 b(up t un! up) — blug, ut, va), 


则 
B2{[e?, e], [vn, an]) 
= v(Vet; Von) — (£2, V un) + (qh, Ve?) 
hk иһ hg 
E У; [x (2 + Vp, Van) + 2. Sv =É (và, иһ), [vOnvn]E) g ) 
KET? EC 
hi. уи" à 
£ > = (= + Vp ,Vgn ) 
KET? 
siu au ajo Y; (m Cy ) — b(un-l al uy) 
Е т" а h Kein 8v т" › V qh K h 3th: ) ЧА 
hA -— 
i De SE (eve Vas). = (u^, u” ‚®һ) + blun uy -i un) + b(ug 1 uR, Uh), 
KeT? 
那么 


Bo(let, eg]; lv, qj) 


Nn $,n—1 
= Bolet, ez] v. 4) — Ba (lets e3], [ns an]) + (“A — диын, vn) 


) HOUR uR vn) — bu t, up, vn) 


u” me 
+ >` SE (T e Vp Vas) и", u”, vn) + bug, up 1 vn) 
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= У { (— vAu” + Vp? + (и" - У)и", о — un) g + (vu — Vpg,v — оһ)к 
KET? 


n 
un 一 


utl 1 B 
(дш = = A e) -3 x (v8 ug]g, v — vn) } + (ип lu uk) 
ECOK 


h2 е u? — yet ss js А 
E ex (E(f = = = Veh, Van) y 7 pa A [идип], [va va]g) g ) 


—b(u^, u^, v) + b(ug, un l un) — b(ug ug un) — (q — qh, V: U&)K 


2 т _ nl 
=- Y {E(P -A ута) - у, PE ((vdnufle,[vOnrnle) в} 


KeTp ECOK 
+ у, (e иһ uh” _ vot, — w) + (f" — f, — uh) 
h n ho K h* K 
KET? 
utt 1 
(аш — h UR v) 7 (7 a5 V uk — 3 > (дли), v — vn) } 


ECOK 


—b(u^, u^, v) + b(ug, ug vx) + b(ug t, uR, vn) 一 b(ug t up lun). 
НХ оһ = Inv, gn = Пд, H Schwarz 不 等 式 与 逆 不 等 式 及 庞 加 莱 不 等 式 知 


llvall < c(h + Diol, 


而 
blug up, un) < eluslilun ilvnli, 
lunli € ch? va] € ch-*(A + Dlvli < eh™ vh, 
所 以 
|b(u, ип оь) € ch ин ив hloh. 
[HJ 2E 
blug, uk, оь) € ch |ushlus hloh, 
Ib(u271, ug, 0) € ch "jug [slug hloh, 
则 


n. n-—1 
Bo (E.B) bd) < ef J [02 + Ма" — fell lee 
KETE 


А 


_ Е 1/2 
+ и" + А2013 (ид) (leh + lal). 
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而 Bolle], єў], [v, q]) 满足 infsup 条 件 (2)， 则 


3 


u? — u^ 
ТЇ + 121° «e JO Ге) + HEE” — feli + ||P — Әл 
KETE 


+ ju” и" + ^"? ( lugi? илт? + jug! I1), 


n (| on|2 nll2 n n2 2 рт nj2 иһ -up ? 
7" (16212 + 1621) < er { > | (nk) +hxllf - feli + |F - Au arm M 
Kern 
+ lu" [u^ B +A? ( lullu? + lu I) }, 
上 述 不 等 式 左右 两 端 关于 n 求 和 ， 即 得 (16) 式 成 立 。 证 毕 


6 ”数值 算 例 


我 们 将 通过 方 腔 流 问题 的 数值 试验 来 验证 本 文 提 出 的 自 适应 变 分 多 尺度 方法 的 有 效 性 。 
$a = 1, аз = 0.5, ñ, = 8 = 0.02， 即 我 们 只 考虑 空间 上 加 密 ， 时 间 步 长 增 大 或 缩小 的 情 
况 。 我 们 取 区 域 Q = (0,1) х (0,1), #77] = 0， 在 区 域 的 上 边界 上 wi = 1, wz = 0， 在 其 它 的 
WHER u = 0, u> = 0. 

结合 自 适应 算法 ， 本 文 给 出 了 Re = 1000 Ft = 50s 时 的 数值 模拟 结果 ， 如 图 1 至 图 2 所 
示 。 从 图 1 我 们 发 现 ， 区 域 中 产生 了 大 涡 ， 同 时 区 域 的 左 、 右 下 方 拐角 处 也 出 现 了 小 涡 。 这 些 
数值 结果 与 文献 [12] 中 的 结果 是 一 致 的 。 图 3 与 图 4 表示 了 Re = 1000 时 两 个 不 同时 刻 的 网 格 
剖 分 情况 ， 从 中 我 们 可 以 看 出 ， 随 着 时 间 的 推移 ， 区 域 顶 端 左 、 右 两 个 抛 角 处 网 格 剖 分 逐渐 加 
密 。 图 5 表示 Re = 1000 时 的 时 间 步 长 变化 情况 ， 体 现 了 自 适 应 算法 的 过 程 。 
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^ 
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图 1: Re = 1000, t= 50s 速度 向 量 图 图 2: Re = 1000, t = 50s 压 力 等 高 线 图 
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图 3: Re = 1000, t= 1s 网 格 训 分 情况 图 4: Re = 1000, t = 50s 网 格 剖 分 情 ; 
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图 5: Re =1000 时 间 步 长 变化 情况 
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Multiscale Adaptivity for the Time-dependent 
Navier-Stokes Equations 


MU Bao-ying, HOU Yan-ren, ZHANG Yun-zhang 


(Faculty of Science, Xi'an Jiaotong University, Xi'an 710049) 


Abstract: In this paper we present a stabilized finite element method for the time-dependent Navier- 
Stokes equations based on the variational multiscale theory. Using a posteriori error estimator which 
is equivalent to the actual error and combining with the adaptive algorithm, we obtain a multiscale 
adaptive method for the time-dependent Navier-Stokes Equations. This method is simple, intuitive, 
and easy to implement. The numerical experiment demonstrates the efficiency of the method. 
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